To appear in the September 21, 2000 issue of Nature 



The Accelerations of Stars Orbiting 

the Milky Way's Central Black Hole 

O" 

O 

O . A. M. Ghez*, M. Morris*, E. E. Becklin*, A. Tanner*, & T. Kremenek* 

(N 

P^; *Dept. of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1562 

(D 



(N 



Recent measurements^ ^, of the velocities of stars near the center of 
the Milky Way have provided the strongest evidence for the presence of a 



K*" \ supermassive black hole in a galaxy^, but the observational uncertainties 

E^ \ poorly constrain many of the properties of the black hole. Determining the 

(^ ' accelerations of stars in their orbits around the center provides much more 

O' 

O ■ precise information about the position and mass of the black hole. Here we 

Q^' report measurements of the accelerations for three stars located ~ 0.005 pc 

^ ■ (projected on the sky) from the central radio source Sagittarius A* (Sgr A*); 

cd . these accelerations are comparable to those experienced by the Earth as it 

■ , . orbits the Sun. These data increase the inferred minimum mass density in the 

H , central region of the Galaxy by an order of magnitude relative to previous 

results and localized the dark mass to within 0.05 ± 0.04 arcsec of the nominal 

position of Sgr A*. In addition, the orbital period of one of the observed stars 

could be as short as 15 years, allowing us the opportunity in the near future to 

observe an entire period. 

In 1995, we initiated a program of high resolution 2.2 ^m (K band) imaging of the 
inner 5 arcsec x 5 arcsec of the Galaxy's central stellar cluster with the W. M. Keck 
10-m telescope on Mauna Kea, Hawaii (1 arcsec = 0.04 pc at the distance to the Galactic 
Center, 8 kpc, ref. 6). From each observation, several thousand short exposure (texp=0.137 
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sec) frames were collected, using the facility near- infrared camera, NIRC^'^, and combined 
to produce a final diffraction-limited image, having an angular resolution of 0.05 arcsec. 
Between 1995 and 1997, images were obtained once a year with the aim of detecting the 
stars' velocities in the plane of the sky. The results from these measurements are detailed 
in ref. 4; in summary, two-dimensional velocities were measured for 90 stars with simple 
linear fits to the positions as a function of time. These velocities, which reach up to 1400 
km/sec, implied the existence of a 2.6 x IO^Mq black hole coincident (±0.1 arcsec) with the 
nominal location of Sgr A* (ref. 9), the unusual radio source^°'^^'^^ long believed to be the 
counterpart of the putative black hole. 

The new observations presented here were obtained several times a year from 1997 to 
1999 to improve the sensitivity to accelerations. With nine independent measurements, 
we now fit the positions of stars as a function of time with second order polynomials (see 
Fig. 1). The resulting velocity uncertainties are reduced by a factor of 3 compared to our 
earlier work, primarily as a result of the increased time baseline, and, in the central square 
arcsecond, by a factor of 6 that presented in ref. 3, due to our higher angular resolution. 
Among the 90 stars in our original proper motion sample^, we have now detected significant 
accelerations for three stars, SO-1, SO-2 and SO-4 (see Table 1); specifically, these are 
the sources for which the reduced chi-squared for a quadratic fit is smaller than that for 
the linear model of their motions by more than 1. These three stars are independently 
distinguished in our sample, being among the fastest moving stars (v = 560 to 1350 km/sec) 
and among the closest to the nominal position of Sgr A* (riggs = 0.004 to 0.013 pc). With 
accelerations of 2-5 milli-arsecond/yr^, or equivalently (3 — 6) x 10~^ km/sec^, they are 
experiencing accelerations similar to the Earth in its orbit about the Sun. 

Acceleration vectors, in principle, are more precise tools than the velocity vectors for 
studying the central mass distribution. Even projected onto the plane of the sky, each 
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acceleration vector should be oriented in the direction of the central mass, assuming a 
spherically symmetric potential. Thus, the intersection of multiple acceleration vectors is 
the location of the dark mass. Figure 2 shows the acceleration vector's direction for the 
three stars. Within 1 a, these vectors do indeed overlap, and furthermore, the intersection 
point lies a mere 0.05 ± 0.03 arcsec East and 0.02 ± 0.03 arcsec South of the nominal 
position of Sgr A*, consistent with the identification of Sgr A* as the carrier of the mass. 
Previously, with statistical treatments of velocities only, the dynamical center was located 
to within ± 0.1 arcsec (la) of Sgr A*'s position^. This velocity-based measurement is 
unaffected by the increased time baseline, as its uncertainty is dominated by the limited 
number of stars at a given radius. Therefore, the accelerations improve the localization of 
our Galaxy's dynamical center by a factor of 3, which is critical for reliably associating any 
near-infrared source with the black hole given the complexity of the region. 

Like the directions, the magnitudes of the acceleration vectors also constrain the 
central black hole's properties. In three dimensions, the acceleration and radius vectors 
{a^-D and r^-n, respectively) provide a direct measure of the enclosed mass simply by 
a^-D = G X M/r^_£,'^. For a central potential, the acceleration and radius vectors are 
co-aligned and, with a projection angle to the plane of the sky 6, the two dimensional 
projections place the following lower limit on the central mass: Mcos^{6) = a2-D x "'"i-z?/^- 
This analysis is independent of the star's orbital parameters, although 6 is in fact a lower 
limit for the orbital inclination angle. Figure 3 shows the minimum mass implied by each 
star's two-dimensional acceleration as a function of projected radius, along with dashed 
curves displaying how the implied mass and radius grow with projection angle. For each 
point, the uncertainty in the position of the dynamical center dominates the minimum 
mass uncertainties. Also plotted are the results from the statistical analysis of the velocity 
vectors measured in the plane of the sky, which imply a dark mass of 2.3 — 3.3 x IO^Mq 
inside a radius of 0.015 pc"^'^. If, as has been assumed, this dark mass is in the form of 
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a single supermassive black hole, the enclosed mass should remain level at smaller radii. 
Projection angles of 51 - 56 and 25 - 37 degrees for SO-1 and SO-2, respectively, would yield 
the mass inferred from velocities and place these stars at a mere ~0.008 pc (solid blue 
portions of the limiting mass curves in Figure 3), thus increasing the dark mass density 
implied by velocities by an order of magnitude to 8 x 10^'^Mq/pc^. With smaller projection 
angles, these two stars also allow for the enclosed mass to decrease at smaller radii, as would 
occur in the presence of an extended distribution of dark matter surrounding a less massive 
black hole^^^^^. In contrast to the agreement between the mass distribution inferred from 
the velocities and the accelerations for SO-1 and SO-2, the minimum mass implied by SO-4's 
acceleration is inconsistent at least at the la level. We note, in support of the validity of 
SO-4's acceleration vector, that its orientation is consistent with the intersection of the SO-1 
and SO-2 acceleration vectors (see Fig. 2). Nonetheless, continued monitoring of SO-4 will 
be important for assessing this possible discrepancy. Overall, the individual magnitudes 
of the three acceleration vectors support the existence of a central black hole with mass 
~ 3 X 10^ Mq. 

With acceleration measurements, it is now possible to constrain the individual orbits. 
We checked the conclusions of the previous studies^''^ by assuming that the stars are bound 
to a central mass of 2.3 — 3.3 x IO^Mq located within 0.03 arcsec of the nominal position of 
Sgr A*. Excellent orbital fits were found for SO-1 and SO-2 for the entire range of masses 
and for true focii within 0.01 arcsec of the nominal position of Sgr A*, suggesting that we 
now have comparable accuracy in determining the IR location of Sgr A* (ref. 9) and the 
dynamical center of our Galaxy (see Fig. 4). The orbital solutions for these stars have 
eccentricities ranging from (circular) to 0.9 for SO-1 and 0.5 to 0.9 for SO-2 and periods 
ranging from 35 - 1,200 and 15 - 550 years, respectively, raising the possibility of seeing 
a star make a complete journey around the center of the Galaxy within the foreseeable 
future. Although the fits are not yet unique, they impose a maximum orbital distance 
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from the black hole, or apoapse, of 0.1 pc. This suggests that SO-1 and SO-2 might have 
formed locally. If these stars are indeed young^, then their small apoapse distance presents 
a challenge to classical star formation theories in light of the strong tidal forces created 
by the central black hole and might require a collisional or compressional star formation 
scenario^^'^^. However, dynamical friction may by able to act on a short enough time scale 
to bring these stars in from a much larger distance^*^. More accurate orbital parameters 
are needed to fully address these problems and others such as the distance to the Galactic 
Center^^. The determinations of these parameters will be considerably improved when 
radial velocities are measured for these stars using adaptive optics techniques, as the current 
solutions based on the proper motion data alone predict radial velocities ranging from 200 
to 2000 km/sec. 



-6- 

REFERENCES 

1. Eckart, A. & Genzel, R. Stellar proper motions in the central 0.1 pc of the Galaxy. Mon. 

Not. R. Astron. Soc, 284, 576-598 (1997) 

2. Genzel, R., Eckart, A., Ott, T. and Eisenhauer, F. On the nature of the dark mass in the 

centre of the Milky Way. Mon. Not. R. Astron. Soc, 291-234, 219 (1997) 

3. Genzel, R., Pichon, C., Eckart, A., Gerhard, O.E., and Ott, T. Stellar Dynamics in the 

Galactic centre: Proper Motions and Anisotropy. Mon. Not. R. Astron. Soc, in 
press (2000) 

4. Ghez, A. M., Klein, B. L., Morris, M., and Becklin, E. E. High Proper-Motion Stars 

in the Vicinity of Sagittarius A*: Evidence for a Supermassive Black Hole at the 
Center of Our Galaxy, Astrophys. J., 509, 678-686 (1998) 

5. Maoz, E. 1998, Dynamical Constraints on Alternatives to Supermassive Black Holes in 

Galactic Nuclei. Astrophys. J., 494, L181-L184 (1998) 

6. Reid, M. J. The distance to the center of the Galaxy. Ann. Rev. Astron. Astrophys. 31, 

345-372 (1993) 

7. Matthews, K. and Soifer, B. T. The Near Infrared Camera on the W. M. Keck Telescope. 

in Astronomy with Infrared Arrays: The Next Generation (ed. I. McLean) 239-246 
(Kluwer Academic Publications, Dordrecht, Astrophysics and Space Sciences Library, 
190, 1994) 

8. Matthews, K., Ghez, A. M., Weinberger, A. J., and Neugebauer, G. The First 

Diffraction-Limited Images from the W. M. Keck Telescope. Proc. Astron. Soc. 
Pac, 108, 615-619 (1996) 

9. Menten, K. M., Reid, M. J., Eckart, A., & Genzel, R. The Position of Sgr A*: Accurate 

Alignment of the Radio and Infrared Reference Frames at the Galactic Center. 



-7- 

Astrophys. J., 475, L111-L114 (1997) 

10. Lo, K.Y., Shell, Z.-Q., Zhao, J.-H., & Ho, P.T. Intrinsic Size of SGR A*: 72 

Schwarzschild Radii. Astrophys. J., 508, L61-L64 (1998) 

11. Reid, M.J., Readhead, A.C.S, Vermeulen, R.C, Treuhaft, R.N. The Proper Motion of 

Sagittarius A*. I. First VLBA Results. Astrophys. J., 524, 816-823 (1999) 

12. Backer, D.C, & Sramek, R.A. Proper Motion of the Compact, Nonthermal Radio 

Source in the Galactic Center, Sagittarius A*. Astrophys. J., 524, 805-815 (1999) 

13. Salati, P. & Silk, J. A Stellar Probe of Dark Matter Annihilation in Galactic Nuclei. 

Astrophys. J., 338, 24-31 (1989) 

14. Tsiklauri, D. & Viollier, R.D. Dark Matter Concentration in the Galactic Center. 

Astrophys. J., 500, 591-595 (1998) 

15. Gondolo, P. & Silk, J. Dark matter annihilation at the Galactic Center. Phys. Rev. 

Lett, 83, 1719-1722 (1999) 

16. Torres, D. P., Capozziello, S., & Lambiase, G. A supermassive boson star at the galactic 

center? Astrophys. J., in press (2000) 

17. Romanowsky, A., and Kochanek, C. to appear in proceedings of "Dynamics of Star 

Clusters and the Milky Way" (2000) 

18. Sanders, R. H. The case against a massive black hole at the Galactic Centre. Nature, 

359, 131-132 (1992) 

19. Morris, M., Chez, A.M., & Becklin, E.E. The galactic center black hole: clues for the 

evolution of black holes in galactic nuclei. Advances in Space Research, 23, 959-968 
(1999) 

20. Gerhard, O. The Galactic Center He I Stars: Remains of a Dissolved Young Cluster? 

Astrophys. J. Lett., submitted (2000) 



21. Salim, S., & Gould, A. Sagittarius A* "Visual Binaries": A Direct Measurement of the 

Galactocentric Distance. Astrophys. J., 523, 633-641 (1999) 

22. Babu, G. J., & Feigelson, E. D. Astrostatistics. Chapman & Hall, London, UK (1996). 

Acknowledgments This work was supported by the National Science Foundation and the 
Packard Foundation. We are grateful to James Larkin for swapping telescope time; Ortwin 
Gerhard, Mike Jura, and Alycia Weinberger for useful input; and telescope observing 
assistances Joel Aycock, Teressa Chelminiak, Gary Puniwai, Chuck Sorenson, Wayne Wack, 
Meg Whittle, and software/instrument specialists Al Conrad and Bob Goodrich for their 
help during the observations. The data presented here were obtained at the W.M. Keck 
Observatory, which is operated as a scientific partnership among the California Institute 
of Technology, the University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous financial support of 
the W.M. Keck Foundation. 

Correspondence and requests for materials should be addressed to A.G. (e-mail: 
ghez@astro.ucla. edu) . 



This manuscript was prepared with the AAS IAT[t;X macros v4.0. 



-9- 



Table 1. Measurements for Stars with Significant Accelerations 

SO-1 (SI) SO-2 (S2) SO-4 (S8) 

Radius from Sgr A*-Magnitudc (milli-pc) 4.42 ± 0.03 5.83 ± 0.04 13.15 ± 0.04 

Radius from Sgr A*-Position Angle (degrees) 290.1 ± 0.7 3.1 ± 0.4 117.3 ± 0.2 

Velocity-Magnitude (km/sec) 1350 ± 40 560 ± 20 990 ± 30 

Velocity-Position Angle (degrees) 168 ±2 241 ± 2 129 ± 2 

Acceleration-Magnitude (milli-arcsec/yr^) 2.4 ± 0.7 5.4 ± 0.3 3.2 ± 0.5 

Acceleration-Position Angle (degrees) 80 ± 15 154 ± 4 294 ± 7 

Note. — 1 - The primary nomenclature adopted here and in the text is from Ghez et al. (1998); in 
parentheses star names from Eckart & Gcnzel (1997) are also given; 2 - all quantities listed are derived 
from the second order polynomial fit of the data and are given for epoch 1995.53; 3 - all position angles 
are measured East of North; 4 - all uncertainties are la and are determined by the jackknife resampling 
mcthod^^; 5 - the radius uncertainties listed include only the relative positional uncertainties and not the 
uncertainty in the origin used (the nominal position of Sgr A*; see text for offset to measured dynamical 
center). 



-10- 

Fig. 1. — East-West [left] and North-South [right] positional offsets from the nominal location 
of Sgr A* vs. time for SO-l[top], SO-2 [middle], & SO-4 [bottom]. The offset range shown is 
scaled to the points in each plot and therefore varies from ~ 0."04 to ~ 0."15. Each of these 
stars is located with a precision of ~ 1 - 5 milli-arcseconds in the individual maps and the 
alignment of these positions, which is carried out by minimizing the net displacement of all 
stars as described in Ghez et al.^, has an uncertainty of ~ 3 milli-arcseconds based on the half 
sample bootstrap resampling method^^ (the alignment uncertainty is at a minimum at the 
center of the map and grows linearly with radius). These two uncertainty terms are added 
in quadrature; the results are used in the fitting process and depicted as errorbars here. In 
each plot, the solid line shows the second order polynomial used to derive the acceleration 
term. These plots demonstrate that we are able to measure, for the first time, accelerations 
of 2 to 5 milli-arcsecond/yr^ (0.3-0.6 cm/s^) for stars orbiting a supermassive black hole. 

Fig. 2. — The acceleration uncertainty cones and their intersections. The cones' edges 
represent the directions for which the accelerations deviate by 1 o" from their best fit values 
and their vertices are the time-averaged positions, measured relative to the nominal position 
of Sgr A*, rather than the positions listed in Table 1. If the accelerations are caused by 
a single supermassive black hole, or even a spherically symmetric mass distribution, these 
vectors should intersect at a common location, the center of the mass. The la intersection 
point lies 0.05 ± 0.03 arcsec East and 0.02 ± 0.03 arcsec South of the nominal position of 
Sgr A*. The existence of an intersection point suggests that there is indeed a common origin 
for the acceleration and pinpoints the location of the black hole to within 0.03 arcsec. 

Fig. 3. — The minimum enclosed mass implied by each star's acceleration measurement vs. 
projected distance from the newly determined dynamical center position. If the projection 
between the plane of the sky and both the radius and the acceleration vectors is 6, then 
the true mass increases as l/cos^(^) and the true acceleration and radius vectors increase 
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as l/cos(^). Also plotted is the mass range implied from a statistical analysis of the ~ 100 
velocity vectors that have been measured in the plane of the sky'^'^. 

Fig. 4. — The measured motion of SO-1, SO-2, & SO-4 and several allowed orbital solutions. 
Only the measurements obtained in ~June of each of the 5 years are shown. SO-1 & SO-2 
are moving clockwise about Sgr A* and SO-4 is traveling radially outward. In the orbital 
analysis, two constraints are applied, a central mass of 2.6 x IO^Mq and a true focus located 
at the position of Sgr A*. Displayed are orbital solutions with periods of 17, 80, 505 years 
for SO-2 and 63, 200, 966 years for SO-1. 
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Figure 1. Ghez et al. 
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Figure 2. Ghez et al. 
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Figure 3. Ghez et al. 
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Figure 4. Ghez et al. 



